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This  report  concerns  the  development  of  methods  for  obtaining  the  optical 
constants  of  anisotropic  crystals  of  the  triclinic  and  monoclinic  systems. 

The  principal  method  used,  classical  disnersion  theory,  is  adapted  to  these 
crystal  systems  by  extending  the  Lorentz  line  parameters  to  include  the  angles 
characterizing  the  individual  resonances,  and  by  replacing  the  dielectric 
constant  by  a  dielectric  tensor.  The  sample  crystals  are  gypsum,  orthoclase 
and  chalcanthite.  The  derived  optical  constants  are  shown  to  be  suitable  for^ 
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20. ^-Abstract 


^jnodelinq  the  optical  properties  of  particulate  media  in  the  infrared  spectral 
reqion.  For  those  materials  where  suitable  size  single  crystals  are  not 
available,  an  extension  of  a  previously  used  method  is  applied  to  alabaster, 


a  polvcrvstall ine  material  of  the  monoclinic  crystal  system. 
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INTRODUCTION 


Computer  modeling  of  the  optical  properties  of  particulate  materials  is  a 
valuable  technique  that  is  useful  in  research  concerning  atmospheric  aerosols 
either  natural  or  man-made  (smoke  screens).  It  is  also  important  in  develop¬ 
ment  of  camouflage  paints  or  coatings  and  has  uses  in  the  description  of 
the  optical  and  thermal  properties  of  materials  which  are  important  in  such 
areas  as  insulation,  solar  collection,  remote  sensing  and  related  topics. 

Such  computer  simulations  have  been  seriously  hampered  by  lack  of  the  optical 
constants  (the  real  and  imaginary  parts  of  the  complex  refractive  index), 
particularly  in  the  infrared  spectral  region,  for  many  materials  that  may 
prove  to  be  of  technical  significance. 

A  considerable  body  of  recent  research  has  gone  into  the  development  of 

1  2  3 

methods  for  extracting  the  optical  constants  of  solids  ’  ’  .  At  the  time 

of  the  inception  of  this  work  the  status  of  the  problem  was  that  isotropic 

materials  such  as  amorphous  solids  and  cubic  crystals  were  well  described 

by  a  complex  refractive  index  which  varies  with  the  frequency  of  the 

radiation.  That  refractive  index  could  be  obtained  by  well  known  techniques 

involving  the  measurement  of  the  reflectance  spectrum  of  such  materials 

4  5 

either  at  two  angles  two  polarizations,  or  a  combination  of  transmittance 
and  reflectance^.  It  would  normally  be  thought  that  in  order  to  obtain  the 
values  of  two  quantities  such  as  n,  the  real  and  k,  the  imaginary  part  of 
the  complex  refractive  index,  m  at  any  frequency  two  measurements  would  be 
needed  as  just  discussed.  However,  use  of  a  single  measurement  (at  each 
frequency)  has  been  employed  with  considerable  success  in  recent  years  (i.e. 
using  classical  dispersion  theory)^.  In  brief,  by  invoking  an  assumed 
spectral  line  shape  (the  Lorentz  shape)  one  can  make  a  single  measurement 
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of  reflectance  or  of  transmittance  over  a  considerable  spectral  band  and 
fit  the  data  to  a  set  of  dispersion  parameters  (Lorentz  line  parameters) 
either  using  a  reflection  or  transmission  method. 

8  9 

Another  method  has  had  much  usage  over  the  years,  i.e.,  Kramers  Kronig  ’ 
analysis.  This  method  involves  a  relationship  between  n  and  k,  but  requires 
that  the  spectrum  be  measured  over  all  frequencies  or  that  an  adequate 
estimate  be  made  of  the  high  and  low  frequency  values.  These  values  at  wave- 
:er.gths  or  frequencies  inaccessible  to  measurement  must  be  known  in  order 
that  the  values  within  the  spectral  range  of  interest  be  not  distorted.  This 
is  not  true  of  the  classical  dispersion  parameter  method,  as  it  has  been 
shown^  that  regardless  of  the  values  outside  the  range  of  interest,  if  the 
reflection  bands  are  well  fit  within  a  given  spectral  range,  the  optical 
constants  obtained  are  appropriate  for  that  range. 

Another  advantage  of  the  classical  dispersion  theory  approach  to  extraction 
of  optical  constants  is  the  convenience  for  computations  of  storing  the 
information  in  a  small  file  of  3N  +  1  parameters  where  N  is  the  required 

number  of  oscillators.  Such  tables  of  e  the  dielectric  constant  at  hi qh 

00 

frequency  v  k  the  frequency,  the  strength  and  y  k  the  damping  for  each 
of  the  independent  oscillators  are  very  small  relative  to  the  multitude  of 
values  of  n  and  k  required  to  make  computations  at  a  large  number  of  spectral 
frequencies. 


The  optical  constants  are  calculated  as  needed  at  any  frequency  by  the  use 


of  the  Lorentz  formula 


N 

2 

k= 


where  e  is  the  comolex  dielectric  constant.  T'/nicallv,  snr»ow|v»t  ^ess  than 
?0  lines,  occasionally  as  few  as  1  or  ?  can  he  used  to  reoresent  the  entire 
spectrum  throunhout  the  infrared  renion. 
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It  has  been  shown  that  in  cases  where  the  absorption  coefficient  is  large 

2  1 1 

the  reflection  method  works  quite  well  ’  and  in  cases  where  the  absorption 

coefficient  is  small  compared  to  the  refractive  index  (organics)  a  trans- 

12 

mission  method  using  a  Brewster  angle  of  incidence  is  appropriate  . 

The  optical  constants  so  derived  have  been  used  together  with  our  theoretical 

1 3 

methods  to  predict  or  simulate  the  spectra  of  particulate  materials  ,  lunar 

14  12  15 

soils  ,  paints  and  aerosols  . 

Extensions  of  these  methods  to  uniaxial  symmetry  crystals  have  been  carried 
out7  and  in  practice  turn  out  to  be  straightforward.  Such  crystals  are 
treated  as  if  they  were,  in  fact,  composed  of  two  independent  materials; 
one  corresponding  to  the  ordinary  ray,  one  corresponding  to  the  extraordinary 
ray.  These  separate  refractive  indices  are  used  to  make  suitable  calculations 
of  cross  sections  for  the  materials  involved  and  then  combined  in  a  2  to 
1  ratio  (Ei  c  :E  l ! c ) ^ 6 ’  17  assuming  that  the  particles  are  randomly 
distributed  in  orientation,  c  refers  to  the  crystallographic  high  symmetry 
axis.  This  technique  has  proven  to  result  in  simulated  spectra  very  closely 
approximating  the  measured  spectra.  The  natural  extension  of  this  method 
to  biaxial  crystals  is  obvious  for  those  of  the  orthorhombic  system.  Using 
a  polarizer  and  appropriate  crystal  orientations  one  may  extract  three  sets 
of  optical  constants11  and  S'/hen  theoretical  computations  for  powders  or 
aerosols  are  made  these  are  used  in  ratio  of  1  to  1  to  1  for  (E  II  a, 

E  II  b  and  E  lie  )  again  assuming  no  preferred  orientation. 
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The  present  report  details  how  the  methods  for  obtaining  the  optical  constants 

have  been  extended  to  the  low  symmetry  monoclinic  and  tri clinic  crystal 

systems.  We  were  aided  in  developing  our  methodology  by  finding  two  papers 

18  19 

dealing  with  the  monoclinic  problem  ’  .  These  papers  suggested  a  method  of 

carrying  out  a  dispersion  theory  analysis  for  monoclinic  crystals  using 

measurements  obtained  from  the  a-c  crystallographic  plane.  We  studied  the  work 

of  these  authors  and  combined  their  approach  with  the  approach  that  we  have 

already  taken  with  respect  to  the  simpler  crystal  systems  .  The  results  of 

this  work  applied  to  gypsum  were  published  in  Ref.  2C.  In  that  paper  we  not 

only  derived  the  optical  constants  for  gypsum,  but  showed  that  these  optical 

constants  were  appropriate  for  simulating  the  emittance  spectra  of  gypsum 

powder  using  our  previously  developed  theory.  We  also  used  the  parameters 

18  19  21 

developed  by  the  Russians  ’  ’  for  spodumene  and  showed  that,  as  with 

gypsum  the  optical  constants  are  indeed  suitable  for  modeling  the  emittance 

of  powders.  Further  monoclinic  work  has  recently  been  submitted  for  publica- 
22 

tion  concerning  the  mineral  orthoclase. 

The  most  difficult  problem,  that  of  triclinic  crystals-,  was  attacked  using 

the  insights  gained  by  the  study  of  the  monoclinic  problem.  The  theoretical 

paper  concerning  the  method  of  carrying  out  the  analysis  was,  in  fact,  the 

23 

first  one  published  during  the  work  on  this  contract.  A  reduction  of  that 

24 

method  to  practice  was  later  published  using  Chalcanthite  (CuSO^'5  i^O). 

Thus  the  body  of  our  work  on  this  contract  has  been  described  in  four  publica¬ 
tions. 


4 


A  Arthur  D.  Little,  Inc. 


However,  some  of  the  details  have  not  yet  been  published.  The  following 

section  will  delineate  the  unpublished  work.  That  work  includes  study  of  the 

monoclinic  substance  alabaster,  which  is  polycrystalline  gypsum.  It  is 

intended  as  an  example  of  the  usefulness  of  carrying  out  classical  dispersioi 

analysis  of  materials  for  which  sufficiently  large  single  crystals  are  not 

25 

available.  We  had  previously  carried  out  such  analysis  for  quartzite  but 
needed  to  show  that  the  same  type  of  analysis  could  be  carried  out  for  the 
low  symmetry  crystal  systems. 

In  addition,  we  show  a  comparison  of  the  measured  experimental  emittance  of 
chancanthite  powder  to  the  predicted  spectrum  using  the  optical  constants  of 
that  triclinic  crystal.  This  was  not  included  in  the  paper  on  chalcanthite 
owing  to  the  simplicity  of  the  spectrum,  which  we  felt  would  not  contribute 
materially  to  the  paper. 

Another  biaxial  crystal  powder,  orthorhombic  enstatite  was  measured  and  its 
spectrum  simulated  using  optical  constants  previously  obtained1^. 

Finally,  the  triclinic  problem  as  discussed  in  Reference  27  proved  quite 
difficult  to  implement  in  such  a  way  that  the  convergence  was  entirely  satis 
factory.  While  the  optical  constants  obtained  proved  useful  and  sensible  we 
believe  that  the  convergence  could  yet  be  improved.  We  therefore  did  not 
previously  describe  the  mathematics  of  the  computer  program  in  detail,  but 
feel  that  it  should  be  presented  here,  if  only  to  give  readers  a  chance  to 
improve  on  our  algorithms. 
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Numerical  Method  for  Triclinic  Crystals 


A.  Preface 

In  the  description  of  the  numerical  methods  used  to  fit  the  theoretical  spectra  to  the  observed 
spectra,  it  is  assumed  that  the  reader  is  conversant  with  references  2  and  23.  In  the  present 
work  we  follow  the  scheme  outlined  in  the  latter  paper  with  certain  exceptions,  the  principal 
change  being  that  we  use  normalized  eigenvectors  in  place  of  the  eigenvectors  (a-j,  b^ ), 

(a2,  b2>  defined  by  Equations  (43)  through  (46)  in  reference  23.  We  also  use  numerical 
subscripts  1,2,3  to  designate  x,  y,  z  axes,  respectively. 

B.  Dielectric  Tensor  of  Triclinic  Crystals 
1.  Discussion  of  Eigenvalue  Problem 

In  order  to  apply  our  usual  method  to  find  the  Lorentz  line  parameters  we  must  calculate 
the  derivatives  of  |r,  { |2  +  |r[  2 12  and  hence  r, ,  and  2  with  respect  to  all  the 
parameters.  As  an  instructive  preliminary,  however,  we  shall  study  the  unusual  eigenvalue 
problem. 


exj  =  Pj!j  ,  j-  1.2 


(1) 


which  plays  a  central  role  in  the  derivation  of  the  reflectance  of  triclinic  crystals.  In  equation 
(1)  we  have  set 


'el  1  <M2  <M3 


e  =  |  e, 2  e22  e 


=  23  I  >the  dielectric  tensor 


(2) 


XJ  " 


el3  e!3  e33. 


klj 


'3  j 


,  the  augmented  eigenvector 


(3) 


fi  - 


*  j 


x 

•> 

.0' 


,  the  eigenvector 


(4) 


and  Pj  is  the  eigenvalue,  which  may  be  found  by  solution  of  the  quadratic  equation 
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-  Pj 

*1  2 

*1  3 

*12 

*22  * 

“  Pj  *2  3 

= 

0 

(5) 

*53 

*2  3 

*33 

We  normalize  the  eigenvector  x.  sc  that 

xfj  +  *ij  =  1  16) 


Let  aj  —  e 1 2  e23  ~  *i  3  (e2 2  —  Pj) 
bj  =  e  1 2  c- 1 3  _  *2  3  (e  1 1  —  Pj ) 


Cj  =(*11-  Pj)  (*22  -  Pj)  -  *12 


(7) 

(8) 
(9) 


Then  x  s  -  = 


fa’  +  b  ?)* 


x2j  = 


bi 


(a2  +  bf)* 


°i 


^3  j 


(a?  +  b2)* 


(10) 

(ID 


(12) 


On  account  of  the  symmetry  of  e,  the  vectors  &  ?2  are  orthogonal.  We  take  the  sign 


/IL  Arthur  D.  Little,  Inc. 


7 


of  the  second  eigenvector  so  that 


Xj  i  ^22*^12“  — ^21  (13), (14) 

This  choice  of  sign  convention  is  arbitrary,  but  it  is  imperative  to  stick  to  one  convention. 

Ours  makes  the  determinant 

Xn  x,2 

=  1  (15) 

X2  1  X2  2 

instead  of  —1. 

In  terms  of  the  components  of  the  normalized  eigenvectors  we  find  the  amplitude  reflectances 
to  be 


where  qj  =  pM2  [ Re(qj)  >  0]  (18) 


Before  we  leave  the  eigenvector  problem,  we  calculate  derivatives  of  the  eigenvalues  and 
augmented  eigenvector  Xj  with  respect  to  the  elements  of  e  in  anticipation  of  our  later 
need  for  them. 
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To  start,  take  the  differential  of  Eq.  (1 ) 


But 


so 


dxTJ  =  0 


?  edx  =  -  Xj  (e13dx,  +  e:3dx,  +  e33dx3) 


(27) 

(28) 


l*1 

b 

i 

C- 


Si' 

V 

■w 

K 


edx  +  dex  =  dpf  +  pdf 


(19) 


& 


The  subscript  j  has  been  omitted,  but  remains  understood.  Multiply  by  fT  (the  superscript  T 
stands  for  "transpose”) 


f1  edx  +  fTdex  =  dpfTf  +  pfTdf 


Note  fTf  =  1  and  (fT+dfT)  (f+df)  =  1  so 


fTdf  =  0 


(20) 


(21) 


!? 


ir 

u 

! 


Thus  we  have 


Since 


it  follows  that 


fTedx  +  fTdex  =  dp 


x  =  f  + 1  0 
x3 


ef  =  ex  -  e  l  0 


(22) 

(23) 


(24) 


V 

.i., 

£ 

J 

B 

rr 


l  * 

.‘J 

<*, 


or 


Therefore 


/6i  3  \ 

e?  =  Pf  “  (  6: 3  x3 

\  3  / 

fTedx  =  dxTef  =  pdxTf  -  dxT  (  e23  j  x3 

\e33/ 


(25) 


(26) 


h\ 

th 

w 


Because  the  third  element  of  f  is  zero  we  see  that 


e 1 3 x ,  +e33.v-2  +  e33x3  =  0 


(29) 


m 


1} 
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and,  therefore,  e,  3dx,  +  e23dx2  +  e33dx3  =  —  X[deI3  -  x2de23  “  x3de33 


(30) 


so  fTedx  =  x1x3de!  3  +  x2x3de23  +  x3  de33  (31) 

Substitute  this  result  into  Eq.  (22)  to  find 

dp  =  fTdex  +  X!  x,;  de,  3  +  x2x3de2  3  +  x3de33  (32) 

Now  all  we  have  to  do  is  to  let  one  dea^  be  non-zero  at  a  time  to  find  the  partial 
derivative  of  p  with  respect  to  the  element  of  the  dielectric  tensor  e^.  The  results  may 
be  summarized  very  simply. 


^eaa 


3pj 

^eap 


=  2xajxpj  (a  *  P) 


or  even  more  succinctly 

3Pj 

de  o  (2  5 X(xi 


(33) 


(34) 


where  5^  is  the  Kronecker  symbol  (unity  when  <x  =  p,  zero  otherwise) 

We  still  must  calculate  the  derivatives  of  x1 , ,  x2,and  x3 ,  with  respect  to  the  elements  eap. 

We  cannot  solve  Eq.  (19)  for  dx, , ,  dx2  ( ,  dx3 ,  because  the  matrix  of  the  coefficients  of  these 
differentials  is  singular  [Eq.  (5)  ]  ,  but  Eq.  (21)  supplies  an  independent  condition.  Since 
f ,  and  f  2  are  orthogonal,  it  follows  that  dx , ,  and  dx,  ,  are  proportional  to  x ,  2  and  x2 2 , 
respectively.  In  fact  we  may  write 


^xl  1  A ap  xl  2  ^ea)3  ’  ^*2 1  ^ aP  x2 2^ea0  »  ^*3 1  “  ^o^x3  2+^cr)3^ea|3  (35) 

We  also  find  ,  when  we  assume  only  one  de^^O,  that  we  can  write 
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with  ua  =  xal  for  a  =  (3 

ua  “  x0i  ’  up  “  x  ai  for  a  =£  (3 

the  otherwise  undefined  components  of  u  being  zero. 

Now  substitute  Equations  (35)  and  (36)  into  Equation  (19)  and  divide  by  de^.  In 
expanded  form  we  see  that 


Pi  0  0 

0  Pi  0 

0  0  0 


1 1 

el  2 

ei3\ 

1  2 

e22 

e2  3  ) 

1  3 

£23 

=  3pi_ 
3ec$ 


Since  e  x2  =  p2  ?2  »  the  first  two  rows  of  (37)  are  seen  to  yield  the  two  equations  for 
Ao0  and  BaP 


(P2-Pi)x,2  A^  +  6,3  Ba<3  =  xl,  2£> - u, 

3eo0 


(P2  Pi  )  X22  \^  +  £23  BQj]=X2|  — 1 —  U2 


Equation  (38)  can  be  solved  by  Cramer's  rule.  For  the  determinant  A  we  find 

A  =  (P2 — Pi )  e33  x3 ! 

since  e;3xI2  -  e13x22  =  ~  e2  3X2 1  —  €| 3X1  1  =  e3  3X3 1 


hence 


Aa/3 


(e2  3X1  [  —  e  1  3 X2  x  )  9P|  ~  623  U|  +  3U3 

_ dJaP _ 

(P2-Pi)e33  x3, 


xnu,  +  X, , u2  - 


= 


Arthur  D.  Little,  Inc. 


11 


All  that  remains  is  to  evaluate  these  expressions  for  pairs  a  &  j3  to  find 


The  derivatives  are  calculated  from  the  equations 


(43) 


d*i  i  =  A^  x,  2 

deap 


dx2 1  _ 
3eo0 


Aa/3 


x3 


9x3i 


=  A^:<3  2  +  BO0 


(44) 


The  derivatives  of  x3  ,  are  not  going  to  be  needed  later,  but  are  included  here  for  completeness 
of  the  exposition. 
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2.  Calculation  of  Dielectric  Tensor  from  Lorentz  Line  Parameters 


Two  sets  of  coordinate  axes  must  be  defined  in  this  aspect  of  the  problem.  One  is  fixed 
in  the  crystal,  the  other  is  fixed  in  the  interferometer  spectrometer. 

The  axes  of  the  crystal  will  be  called  x,  y  &  z.  Let  the  x-axis  be  identical  with  the  a-axis  of  the 
crystal.  Let  the  y-axis  be  in  the  a-b  plane  but  perpendicular  to  the  x-axis.  The  z-axis  is 
placed  perpendicular  to  both  x-axis  &  y-axis  to  form  a  right  handed  system  (xyz). 

The  axes  in  the  interferometer  are  such  that  the  plane  of  polarization  of  the  light  incident 
on  the  crystal  face  is  along  the  1-axis;  the  light  propagates  along  the  3-axis.  The  2-axis 
is  chosen  perpendicular  to  the  1-  &  3-axes  to  form  a  right  handed  system  (1,2,  3). 

To  model  the  experiment  on  the  crystal,  the  crystal  properties,  expressed  in  the  (x,  y,  z) 
axes  must  r;e  rotated  into  the  (1,  2,  3)  system.  To  do  this  we  must  set  up  a  rotation 
matrix  S  (R  is  reserved  for  reflectance)  with  elements  S  j  chosen  so  that 
/Sn\ 

I  S2 1  I  is  the  unit  vector x  expressed  in  the  interferometer  or  lab  system  (45a) 

vw 

/s.A  /s„\ 

I  S22  land  I  S23  j  are  the  unit  vectors  y  and  z  expressed  similarly  (45b) 

\S3  2/  \S3  3/ 

Now  suppose  we  know  the  dieletric  tensor,  7,  referred  to  the  crystal  axes.  Then  the  tensor 
in  the  lab  system  is 


e  =  Se  ST 

(46) 

06 

and  - =  S 

d€  sT 

(47) 

0p?k 

under  the  assumption  that  the  Lorentz  lines  occur  independently  of  each  other  (instead  of 
in  orthogonal  triads)  we  have,  in  the  crystal  system 
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e  =  e"  + 


VkTLk 


where  e~  has  components  e~x>  e^y  etc.  &  is  a  symmetrical  tensor  representing 

the  dielectric  constant  at  high  frequencies. 


sin  0k  cos  \ 

sin  0k  sin  4>k  1 ,  giving  the  orientation  of  the  dipole 


1  +i\  -  V 

"k  \VJ 


,  the  usual  resonance  formula. 


We  further  define  Lorentz  line  parameters  pgk  as  follows: 


P01  P02  P0  3 


I  P0  2  P04  P05 


P0  3  P05  P06 


Plk  =  Sk 

P2k  =  V'k 


P3  k  =  "k 
P4k  =  0k 


P5k  =  \ 


In  order  to  perform  “Newton’s  Descent”  to  calculate  the  best  parameters  to  minimize  the 
evil20,  we  need  the  matrix  UUT  whose  elements  are 


E3R  (v  p)  3R 

—  {  —  (t'cP) 

9pok  9pXk 
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where  E  is  the  total  number  of  experimental  reflectance  measurements  made  and  2k,  Xx  run 
over  all  possible  values  that  have  been  defined. 


As  before 


3R  = 
aP2k 


3rn  +  ri2 
3P£k 


(54) 


and 


0rij  V  3rij  deaB 

—L  =  2-j  — -  —  (6  terms  :  a/M  1,12,13,22,23,33) 

3pgk  “0  3ea0 


(55) 


Since  rt  L 


(56) 


&  r12  =  2  x, ,  x01  - - 

■  y+q,  l+q2 

we  write 


dJlL: 

Jl!L 

^u  + 

.Oil 

3x2t  3rij 

"  *t"  * 

3q. 

_ L.  + 

dri} 

3q, 

3ea0 

3xu 

dea0 

3x,, 

3ea0  3cM 

dea0 

3q, 

3ea0 

From  (56)  &  (57)  we  see 


QJ 

-1 

It 

4  xi  i 

3ri2 

3xu 

1  +qt 

3Xu 

3rn 

4  *21 

3rl2  _ 

3x2l 

1  +q2 

3x2l 

3rll  . 

2xu 

3r!2  _ 

3q  i 

(1  +  q,)2 

3q, 

-  xi 1  X12 

U  +q,)2 


(57) 


(5[ 


(59) 


(60) 


(61) 
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2  X  21 


^^12  2  Xj 2  ^2j 


3q2  0+q2)2 


3q2  (l+q2)2 


We  have  already  demonstrated  that 


d  xi  l  Aaj3  xi  2 


d  xi  2  Aa0  X2  2 


Aj  j 


2  ■> 
X2  1  -XIl 

P  —  P 
F2  M 


A,  3  —  (2e13  Aj  j  +c23  Ai2)/e3  3 


Ao  t  A. 


A23  ^“^23  A22  "^^13  A12^^33 


A3  3  Xj  j  Xj  2/(P2  —  Pj  ) 


=  '/;  X" 


=  *Ct 


So  the  computation  of  — - —  is  defined.  We  must  therefore  calculate  the  derivatives 

dea(S 


.  For  the  entire  matrix 
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and  the  have  been  found  to  be 


_  xi  1  X2 1 


P2-Pi 


Obviously,  from  (33),  we  have 


(67) 


=  S  ST 

9Pfik  aP2k 


So  we  confine  ourselves  to  the  calculation  of  9T  and  suppose  the  results  will  be  rotated 

8Pfik 


numerically  into  the  (1,2,3)  (lab)  system. 


To  start  with. 


000  I  ’  ~ 


0  1  0\ 

1  0  0  J  etc. 

0  0  0/ 


The  first  three  derivatives  of  the  Lorent2  line  part  of  the  affair  are  unchanged  from  before: 


<vkvk  Lk)  =  Vk  —  *  i:z  1.2,3 


9Lk  1 

—  =  —  Dk  =  1  +,P2ki;  -  W 
aPik  Dk 


9Lk  _  _  Lk 


_  ,  L* 

-  -  z  M  -  -- 


But  -■  (VjjVjJ  Lj.)  =  — 

dp?k  \3pgk 


9  t  /  3vk  r  3vk  \ 

7 - (VkM  =  —  vk  +  vk  -  Lk'  2  =  4-5 
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s 


9v£  =  (cos  P4k  cos  p5k  ,  cos  P4k  sin  p5k  ,  -  sin  p4k) 


=  (-sin  p4k  sin  p5k  ,  sin  P4k  cos  p5k  ,  0  ) 


Obviously  Ov, 


H  V 


i  vT  = 
vk 


A  Fortran  program  embodying  these  methods  was  written  and  used  to  fit  the  observed 
spectra.  Convergence,  however  was  not  satisfactory  for  reasons  that  are  still  unknown. 


All  of  the  algebra  and  the  coding  was  checked  by  independent  computation.  To  perform 
the  checking,  a  program  was  written  for  a  HP41CX  to  calculate  the  reflectance  R  as  a 
function  of  the  Lorentz  line  parameters.  Derivatives  were  checked  by  comparing  values 
obtained  by  the  Fortran  program  with  those  obtained  from  the  HP4ICX  by  numerical 
differentiation. 


/k.  Arthur  D.  Little,  Inc. 


III.  Simulations  of  Powder  Spectra  using  Derived  Optical  Constants 

As  mentioned  above,  the  usefulness  of  our  derived  optical  constants  can  be 
evaluated  by  examining  the  spectra  simulated  with  them  in  comparison  to 
measured  emittance  spectra  of  powders.  Such  comparisons  were  carried  out  in 
References  20  and  22  during  the  contract.  Figure  1  shows  such  a  comparison 
tor  chalcanthite  powder.  This  comparison  shows  good  agreement  between  the 
exper  mental  data  and  the  simulated  spectrum  using  the  three  sets  of  optical 
constants  derived  in  Reference  24.  The  spectrum  is,  however,  so  simple  that 
we  did  not  believe  it  provided  a  convincing  proof  of  the  validity  of  our 
method. 

A  better  case  is  shown  for  enstatite  in  Figure  2.  The  three  sets  of  optical 
constants  for  this  material  had  been  derived  previously  by  classical 
dispersion  theory11.  As  enstatite  is  orthorhombic  it  was  used  as  a  test  of 
our  use  of  a  1:1:1  ratio  of  the  three  sets  of  optical  constants  for  randomly 
oriented  biaxial  powders  where  the  crystal  axes  are  orthogonal. 

Finally,  there  remains  the  problem  of  obtaining  useful  optical  constants  for 

modeling  purposes  with  those  crystals  for  which  a  sufficiently  large  single 

crystal  cannot  be  obtained.  The  method  we  have  chosen  to  attempt  for  such 

materials  consists  of  measuring  the  reflectance  spectrum  of  a  polycrystalline 

sample  and  obtaining  optical  constants  as  if  this  material  was  isotropic.  The 

resulting  values  are  clearly  not  rigcrour  representing  the  surface  (Fresnel) 

averaged  spectrum,  as  this  average  is  differtnt  from  the  average  (cross 

sections)  that  would  be  appropriate  for  modeling  of  particulate  materials. 

2  25 

Nonetheless,  this  idea  has  been  explored  and  found  useful  in  the  past.  . 
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FIGURE  1  EMITTANCE  SPECTRUM  OF  CHALCANTHITE 
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FIGURE  2  EMITTANCE  SPECTRUM  OF  ENSTATITE 
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In  order  to  carry  out  the  same  kind  of  test  of  "pseudo"  optical  constants 
for  a  low  symmetry  crystal  we  chose  to  use  alabaster  which  is  a  poly¬ 
crystalline  form  of  gypsum. 

A  sample  obtained  from  Wards  Natural  Science  Establishment  was  cut  and 

polished  and  its  reflectance  spectrum  obtained  using  our  Digilab  FTS  15C 

20 

Interferometer  Spectrometer  with  the  Harrick  MIRA  attachment.  Rapid 
convergence  to  the  Lorentz  line  parameters  was  obtained.  The  values  for 
alabaster  are  given  in  Table  1.  The.  resulting  fit  is  shown  in  Figure  3. 

In  Figure  4  we  show  the  simulated  spectra  of  the  gypsum  powder  of  Reference 
20  using  both  the  true  optical  constants  derived  from  single  crystal  gypsum 
and  "pseudo"  optical  constants  obtained  from  alabaster.  We  have  used  an 
option  in  our  powder  theory  program  that  allows  for  the  contributions  of 
asperities  on  the  gypsum  particles  here  which  results  in  a  sTghtly  differ¬ 
ent  intensity  than  shown  in  the  theoretical  spectrum  of  Reference  20.  As  is 
evident  the  pseudo  optical  constants  do  not  give  as  good  a  fit  to  powder 
data  as  the  true  optical  constants.  The  nature  of  chat  defect  is  principally 
that  the  simulated  emittance  features  are  much  weaker  than  those  observed. 
This  is  a  similar  result  to  that  obtained  in  our  previous  quartzite  experi¬ 
ment  but  more  striking.  We  believe  the  discrepancy  to  result  from  the 
difficulty  of  obtaining  a  suitably  high  reflectance  on  alabaster  which  is 
hard  to  polish  well  owing  to  its  soft  character.  A  further  discrepancy  is 
the  relative  large  size  of  the  feature  near  600  cm  ^  compared  with  either 
the  experiment  with  gypsum  powder  or  the  simulation  with  the  true  optical 
constants. 
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FIGURE  4  EMITTANCE  SPECTRUM  OF  GYPSUM 


IV. 


Summary 


Methods  for  derivation  of  the  optical  constants  of  materials  belonging  to 
the  triclinic  and  monoclinic  crystal  systems  have  been  developed.  These 
optical  constants  were  shown  to  be  useful  in  modeling  the  emittance  spectra 
of  powders  in  the  infrared  region.  The  rigorous  method,  which  can  only  be 
used  for  single  crystals,  provides  better  results  than  an  approximate  method 
which  must  be  used  if  sufficiently  large  single  crystals  are  not  available. 

The  bulk  of  this  work  has  been  published  or  submitted  for  publication  in 
four  papers.  They  are: 

1.  A.  G.  Emslie  and  J.  R.  Aronson,  Determination  of  the  Complex  Dielectric 
Tensor  of  Triclinic  Crystals:  Theory,  J.  Opt.  Soc.  Am.  73,  916  (1983) 

2.  James  R.  Aronson,  Alfred  G.  Emslie,  Ellen  V.  Miseo,  Emmett  M.  Smith  and 
Peter  F.  Strong,  Optical  Constants  of  Monoclinic  Anisotropic  Crystals: 
Gypsum,  Appl  Opt.  22,  4093  (1983) 

3.  James  R.  Aronson,  Alfred  G.  Emslie  and  Peter  F.  Strong,  Optical 
Constants  of  Triclinic  Anisotropic  Crystals:  Blue  Vitriol, 

Appl.  Opt.  24,  1200  (1985) 

4.  James  R.  Aronson,  Optical  Constants  of  Monoclinic  Anisotropic  Crystals: 
Orthoclase,  submitted  to  Spectrochimica  Acta  (1985) 
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